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Abstract. Application of a diagnostic system to a Helicopter gearbox is presented. THe diagnostic system is a 
nonparametric pattern classifier tHat uses a multi-valued influence matrix (MVIM) as its diagnostic model and 
benefits from a fast learning algorithm that enables it to estimate its diagnostic model from a small number 
of measurement-fault data. To test this diagnostic system, vibration measurements were collected from a 
helicopter gearbox test stand during accelerated fatigue tests and at various fault instances. The diagnostic 
results indicate that the MVIM system can accurately detect and diagnose various gearbox faults so long as 
they are included in training. 


1. INTRODUCTION 


2. THE MVIM METHOD 


Helicopter drive trains are significant contributors to both 
maintenance cost and fright safety incidents. Drive trains 
comprise almost 30% of maintenance costs and 16% of me- 
chanically related malfunctions that often result in the loss 
of aircraft (Chin, 1993). As such, it is crucial that faults 
be diagnosed in-flight so as to prevent loss of lives. 

Fault diagnosis of helicopter power trains is based primar- 
ily on vibration monitoring. As such, considerable effort 
has been directed toward the identification of features of 
vibration that are affected by specific faults (e.g., Pratt, 
1986), and the development of signal processing techniques 
that can quantify such features. The main problem with 
this approach, however, is that due to the complexity of 
helicopter gearboxes and the interaction between their var- 
ious components, the individual vibration features do not 
provide a reliable basis for diagnosis. 

As an alternative to single-feature based diagnosis, fault 
diagnosis based on several features can be performed using 
pattern classification (Pau, 1981). Among the various pat- 
tern classifiers used for diagnosis, artificial neural nets are 
the most notable due to their nonparametric nature (inde- 
pendence of the probabilistic structure of the system) and 
their ability to generate complex decision regions. How- 
ever, neural nets generally require extensive training to 
develop the decision regions (diagnostic model). In cases 
such as helicopter power trains, where adequate data may 
not be available for training, artificial neural nets may mis- 
diagnose the fault. 


Measurements are processed in the MVIM method as il- 
lustrated in Fig. 1. They are usually filtered first to ob- 
tain a vector of processed measurements P, then they are 
converted to binary numbers through a flagging operation 
(i.e., abnormal measurements characterized by ‘1’ and nor- 
mal ones by ‘0’), and finally they are analyzed through 
the diagnostic model. In the MVIM method, flagging is 
performed by a Flagging Unit that is tuned according to 
measures of diagnosability and fault signature variability 
obtained from MVIM so as to improve the fault signa- 
tures. The MVIM method is explained in detail in (Danai 
and Chin, 1991) and (Chin and Danai, 1992). We will only 
discuss the overall concept here for completeness. 
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Fig. 1: Processing of measurements in the MVIM 
method 


2.1. Fault Signature Representation 

Fault signatures in the MVIM method are represented by 
the n unit-length columns V ; 6 'R m of a multi-valued in- 
fluence matrix (MVIM) A: 


In this paper we demonstrate the application of a diagnos- 
tic method that can establish the fault signatures based on 
a small number of measurement-fault data. This method 
utilizes a multi-valued influence matrix (MVIM) as its di- 
agnostic model which provides indices for diagnosability of 
the system and variability of the fault signatures (Danai 
and Chin, 1991). These indices are used as feedback to 
improve fault signatures through adaptation (Chin and 
Danai. 1993). 


A= [V,, V,-, ... V n ] (1) 

where m denotes the number of measurements, and n rep- 
resents the number of faults. Based upon this influence 
matrix, the faults can be ranked according to their pos- 
sibility of occurrence by the closeness of their influence 
vectors V, to the vector of flagged measurements Y (see 
Fig. 2). In the MVIM method, the vector of diagnostic 
certainty measures , X, which ranks the faults according to 
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possibility of occurrence is defined as 

X = {xi, Xjy x»} T = cos {on, Oj, o„} T (2) 

where the Xj represent the individual diagnostic certainty 
measures, and the denote the individual angles between 
the influence vectors V ; and the flagged measurement vec- 
tor Y (see Fig. 2). 



Fig. 2: Schematic of diagnostic reasoning in the MVIM 
method, illustrated in three dimensional space 


The influence vectors in Eq. (1) are not known a-priori and 
need to be estimated. In the MVTM method, the error in 
dia gnosis is used as the basis to estimate/update the in- 
fluence vectors. For this purpose, the fault signatures are 
updated recursively after the occurrence of each fault to 
minimize the sum of the squared diagnostic error associ- 
ated with that fault (Danai and Chin, 1991). 

2.3. Fault Signature Evaluation 

One of the unique features of the MVIM method is its 
ability to evaluate quantitatively the uniqueness and vari- 
ability of fault signatures, so that these quantitative mea- 
sures can be used to improve the flagging operation. In 
the MVIM method, the uniqueness of fault signatures is 
represented by the closeness of pairs of influence vectors, 
and the index of diagnosability is defined to characterize 
the closest pair of fault signatures. 

In the MVIM method, the variability of fault signatures 
is defined by their variance. For this purpose, the vari- 
ance matrix of A is estimated to provide a^ measure of the 
variations of individual components of A. Since in the 
MVIM method the components of A are adjusted recur- 
sively, the variance matrix can be readily estimated during 
training (Chin and Danai, 1993). The index of fault signa- 
ture variability is defined as the largest component of the 
variance matrix, representing the largest variability in the 
components of A. 

2.4. Flagging 

The influence matrix A is estimated based on the values of 
the flagged measurement vector Y. Thus, before the influ- 
ence matrix is used for diagnostic reasoning, the integrity 
of the flagging operation needs to be ensured. Flagging in 
the MVTM method is performed by a Flagging Unit that 
is tuned to improve the diagnosability of the system and 
reduce the variability of the fault signatures (Chin and 
Danai, 1993). 

The Flagging Unit uses a sample set of measurement -fault 
vectors to tune its parameters iteratively. After each pass 
through the training batch, the Flagging Unit counts the 


total number of false alarms and undetected faults, and 
estimates the uniqueness and variability of the estimated 
fault signatures from the current values of the influence 
matrix, so that it can use these measures as feedback in 
the next adaptation round. Adaptation stops when the 
total number of false alarms and undetected faults are 
minimized, and the uniqueness and consistency of fault 
signatures are enhanced. 

The Flagging Unit processes the residuals as follows (see 
Fig. 3). The residual vector P G is first passed through 
Hard-Limiter I (consisting of a vector of m thresholds, hi,, 
i = \ f ... ,m) to produce a binary vector Z G B ™ . This 
vector is then multiplied sequentially by the normalized 
columns of a Quantization Matrix and then thresholded 
by Hard-Limiter II to produce the individual components 
of the flagged measurement vector Y G B™ . Training of 
the Flagging Unit comprises of adjusting the thresholds of 
Hard- Limiters I and II and the Quantization Matrix, as 
explained in detail in (Chin and Danai, 1992). 



Fig. 3: Schematic of the Flagging Unit 


3. EXPERIMENTAL 

Vibration data were collected at NASA Lewis Research 
Center. Various component failures in an OH- 58 A main 
rotor transmission were produced during the experiments 
(Lewicki et al , 1992). The configuration of the trans- 
mission which was tested in the NASA 500-hp Helicopter 
Transmission Test Stand is shown in Fig. 4. The vibration 
signals were recorded from eight piezoelectric accelerome- 
ters (frequency range of up to 10 KHz) using an FM tape 
recorder. The signals were recorded once every hour for 
about one to two minutes per recording (at the tape speed 
of 30 in /sec, providing a bandwidth of 20 KHz). 

A total of five tests were performed, where each test was 
run between nine to fifteen days for approximately four 
to eight hours a day. New components were used at the 
start of ***rifr test. When a component fault was detected 
during a test, it was replaced with a new one for the re- 
mainder of the test. Eleven failures occurred during these 
tests. The status of various faults during the five tests are 
shown in Table 1 where no- fault cases are denoted by Jo- 
in Test $1, spiral bevel pinion failure (J 4 ) is estimated to 
have been present on days 5-9 with sun gear failure (x,) 
also occurring on day 9. No failures occurred in Test #2, 
so all the nine days for this test are marked as normal. 
In Test #3, planet bearing failure (x 2 ) was established 
to have been present on days 3-4 and 11-12, with hous- 
ing crack (r 3 ) and mast bearing failure (x 6 ) occurring on 
day 9 and day 13, respectively. Test #4 is estimated to 
have contained planet bearing failure on days 11-12 and 
sun gear failure on days 14-15. In Test #5, housing crack 
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Fig. 4: Configuration of the OH-58A main rotor trans- 
mission 

is assumed to have been present on days 7-11, with sun 
gear failure (xj) and planet gear failure (x 5 ) occurring on 
days 10-11 and day 11, respectively. 


eters were estimated for each segment and averaged over 
these segments. The data records as well as the param- 
eters obtained from the above processing modules were 
then transferred to a personal computer for further analy- 
sis (Chin, 1992). 

5. IMPLEMENTATION AND RESULTS 

The configuration of the MVIM system as applied to fault 
diagnosis of the OH-58A main rotor transmission is illus- 
trated in Fig. 5. As shown in this figure, two MVIMs 
are trained for each accelerometer. One MVIM to per- 
form detection (i.e., to determine whether a fault has oc- 
curred or not), and a diagnostic MVIM to isolate the fault. 
The 54 parameters obtained from the signal analyzer were 
used to train and test the MVIM system. The detection 
MVIM contained only two columns to characterize the no- 
fault and fault signatures, whereas the diagnostic MVIM 
contained seven columns, one characterizing the no-fault 
signature and the other six representing the signatures of 
individual faults. Note that the two MVIMs can be per- 
ceived as filters with different resolutions. In order to in- 
tegrate the results from the MVIMs associated with the 
eight accelerometers, a voting scheme was utilized. Test 


Table 1 : Association of data from each day of the five tests 
with the no-fault case x 0 and various fault cases. 
The six faults are represented as x\: sun gear 
failure , x 2 ; planet bearing failure , x 3 : housing 
crack, x 4 : spiral bevel pinion failure, x 5 : planet 
gear failure, and x € : mast bearing failure 


Fault Status 

Test #1 Test #2 Test #3 Test #4 Test #5~ 


1 

x 0 

*0 

Xo 

Xo 

Xo 

2 

x 0 

x 0 

Xo 

Xo 

Xo 

3 

x 0 

x 0 

X 2 

Xo 

x 0 

4 

x 0 

Xo 

X2 

Xo 

Xo 

5 

x 4 

x 0 

x 0 

Xo 

Xo 

6 

x< 

Xo 

Xo 

Xo 

Xo 

7 

x 4 

Xo 

Xo 

Xo 

x 3 

8 

x 4 

x 0 

Xo 

Xo 

X 3 

9 

X 4 ,X! 

Xo 

X 3 

Xo 

X 3 

10 



Xo 

Xo 

X 3 ,X] 

11 



x 2 

X 2 

X 3 ,Xi,X 5 

12 



x 2 

x 2 


13 



X$ 

Xo 


14 




x x 


15 




Xl 



4. SIGNAL PROCESSING 

In order to identify the effect of faults on the vibration 
data, the vibration signals obtained from the five tests were 
digitized and processed by a commercially available signal 
analyzer (Stewart Hughes, 1987) with four processing mod- 
ules: (1) Statistical Analysis (STAT), (2) Baseband Power 
Spectrum Analysis (BBPS), (3) Bearing Analysis (BRGA), 
and (4) Signal Averaging Analysis (SGAV). 

For analysis purposes, only one data record per day was 
used for each test. The data records were taken at the 
beginning of the day unless a fault was reported. When a 
fault was detected, the record was taken right before the 
fault incident to ensure that the data record reflected the 
fault. Also, in order to reduce estimation errors, each data 
record was partitioned into sixteen segments, and param- 



A 

X 

Fig. 5: Configuration of the MVIM system as applied to 
the OH-58A main rotor transmission. 


#3 and #4 contained most of the failure modes (i.e., 4 
out of 6). Therefore, the parameters from these two tests 
were used to train the MVIMs. Note that not all of the 
failure modes were included in training, so the test results 
were not expected to be perfect. For training the detection 
MVIMs, only the 19 parameters from the STAT , BBPS , 
and BRGA modules were used. Previous studies on this 
data show that these 19 parameters are adequate for de- 
tection (Chin, 1993). For training the diagnostic MVIMs, 
all of the 54 parameters were utilized. 

The initial values of the detection MVIMs (19 x 2) and di- 
agnostic MVIMs (54 x 7) were each set to 0 (i.e., matrices 
with all zero entries), and the initial values of the Quan- 
tization Matrices (19 x 19 for detection) and (54 x 54 for 
diagnosis) were set to identity matrices. The initial thresh- 
old levels for Hard-Limiter I were set at the mean plus one 
standard deviation of the corresponding parameter, and 
for Hard-Limiter II, they were set at 0.5. The maximum 
number of epochs for training the detection and diagnos- 
tic MVIMs was set to 50. After each epoch, the detec- 
tion/diagnostic performance of MVIMs within the training 
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set was tested. Training was stopped once perfect detec- 
tion/diagnosis was achieved, to avoid overtraining (Hertz 
et aL , 1991). The number of epochs used for individual 
detection MVIMs were: 8, 5, 50, 37, 50, 15, 50, and 50 
for accelerometers #1 to #8, respectively, whereas for di- 
agnostic MVIMs they were: 50, 1, 2, 2, 26, 50, 50, and 
50. Accor ding to the number of epochs used for individual 
MVIMs, it is clear that the detection MVIMs associated 
with accelerometers #3, #5, #7, and #8 did not achieve 
perfect detection within the training set. Similarly, the di- 
agnostic MVIMs associated with accelerometers #1, #6, 
#7, and #8 did not achieve perfect diagnosis within the 
training set. 

The MVIMs trained on Tests #3 and #4 were evaluated 
for all of the tests. For this purpose, the 19 parameters 
from each of the eight accelerometers were first passed 
through the corresponding detection MVIM for all of the 
five tests to reflect the occurrence of faults. Once a fault 
was posted by a detection MVIM, the set of 54 parame- 
ters from that accelerometer was passed through the cor- 
responding diagnostic MVTM. At the final stage, the diag- 
nostic certainty measures obtained from the eight diagnos- 
tic MVIMs were consolidated by the voting scheme, which 
utilized weights reflecting the speed of convergence of in- 
dividual fault signatures during training (Chin, 1993). As 
such those influence vectors which converged faster were 
assigned larger weights and vice versa. Zero weights were 
assigned to the influence vectors which did not converge 
during training, and unity weights to those which con- 
verged within one epoch. 

The diagnostic results obtained from the above diagnostic 
system are shown in Table 2. The results indicate that 
the MVTM system was able to produce perfect diagnostics 
for Tests #3 and #4, on which it was trained, and that it 
diagnosed 88% of the faults in all of the tests. Specifically, 
the results in Table 2 indicate that the MVIM system pro- 
duced two false alarms (on day 4 of Test #1 and day 6 of 
Test #5), and five misdiagnoses (on days 5-8 of Test #1 
and day 11 of Test #5). In addition, this system produced 
equal diagnostic certainty measures for the no-fault case 
and s un gear failure on day 10 of Test #5 and could only 
diagnose one of the faults on day 9 of Test #1, and on 
days 10 and 11 of Test #5. However, it should be noted 
that faults x 4 and i 5 were not included in training, so 
no fault signatures were estimated for them. The correct 
diagnostic rate of MVIM, with these two faults excluded 
would be over 95%, which is quite noteworthy considering 
that the MVTM system was trained on a small number of 
measurement -fault data with very few repetitions of each 
fault. 

6, CONCLUSION 

A fault diagnostic system based on the MVIM method is 
applied to a helicopter gearbox for which only a small num- 
ber of measurement-fault data was available. This diag- 
nostic system utilizes two levels of isolation and integrates 
the results obtained from various accelerometers through 
a voting sch em e. The diagnostic results indicate that the 
MVIM system correctly detected all of the fault incidents, 
with only two false alarms posted. The results further in- 
dicate that this system correctly diagnosed all of the faults 
it was trained for. 


Table 2: Estimated fault status from each day of the five 
tests. Same notations are adopted as in Table 1 


Estimated Fault Status 

Day 

Test #1 

Test #2 

Test #3 

Test #4 

Test #5 

1 

x 0 

xo 

x 0 

x 0 

xo 

2 

" x 0 

x 0 

x 0 

x 0 

Xo 

3 

x 0 

X 0 

x 2 

Xo 

Xo 

4 

x 3 

XQ 

X 2 

Xo 

x 0 

5 

*3 

Xq 

Xo 

Xo 

x 0 

6 

*3 

Xo 

Xo 

Xo 

*6 

7 

*3 

Xo 

Xo 

Xq 

*3 

8 

X 3 

x 0 

Xo 

Xq 

X 3 

9 

*1 

Xo 

X$ 

X 0 

Xz 

10 



Xo 

Xo 

Xo,x x 

11 



x 2 

X 2 

Xj, x 6 

12 



X 2 

X2 


13 



Xo 

xo 


14 




X\ 


15 




X\ 
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